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Abstract:  This  study  examines  the  moisture- 
tension  constituent  relationship  of  ice  near  0°C  by 
measuring  its  moisture  characteristic  curve.  In  this 
novel  approach,  we  held  ice  samples  in  a  Tempe 
cell  submerged  in  a  cold  bath  and  removed  water 
from  the  ice  veins  through  incremental  increases 
in  air  pressure.  Observed  air  entry  and  water  re¬ 
entry  pressures  generally  support  an  idealized  flow 


path  model  of  alternating  vein  segments  and 
nodes.  Re-entry  of  water  into  the  air-filled  ice  sam¬ 
ples  occurred  slowly  over  a  period  of  weeks,  sug¬ 
gesting  low  permeabilities  on  the  order  of  1  x 
m2.  We  conclude  that  ice  desaturates  and  rewets 
in  a  manner  similar  to  soils,  with  entry  pressures 
on  the  order  of  silts. 


Cover:  Ice  sample  in  a  Tempe  cell  after  being  rotted  for  several  months  in  a 
refrigerated  bath  held  near  0°C. 
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Capillary  Tension  in  Rotting  Ice  Layers 

RACHEL  E.  JORDAN  AND  JEFFREY  A.  STARK 


INTRODUCTION 

Ice  layers  within  snow  or  at  the  ground  interface 
impede  and  divert  the  flow  of  water,  thus  leading  to 
ponding  and  flooding  during  rain-on- snow  events  or 
during  springtime  thaws.  Thawing  of  the  basal  ice  layer 
affects  terrain  and  road  mobility.  We  must  therefore 
question  how  ice  behaves  in  contact  with  soil  and  in 
the  context  of  layered  porous  media.  An  upgraded  ver¬ 
sion  of  the  SNTHERM  snow  model  (Jordan  1991,  Jor¬ 
dan  et  al.  1999)  describes  movement  of  water  under 
thermal  and  liquid  gradients  in  frozen  and  unfrozen  soil. 
In  this  model  (SLTHERM),  temperature,  pressure,  ice, 
and  liquid  are  coupled  through  constituent  relationships 
that  are  universally  used  for  soil  (e.g..  Black  and  Tice 
1989,  Black  1990).  While  water  flow  through  glaciers 
is  often  treated  as  Darcian  saturated  flow  (e.g.,  Colbeck 
1976,  Lliboutry  1996),  here  we  consider  whether  con¬ 
stituent  relationships  similar  to  soil  can  be  applied  to 
ice.  We  ask  the  fundamental  question  whether  air  and 
liquid  pathways  can  co-exist  in  ice,  and  if  so,  what 
configurations  these  phases  take.  In  this  paper  we  con¬ 
centrate  on  the  moisture-tension  constituent  rela¬ 
tionship,  which  characterizes  the  air-liquid  proportion 
in  the  pore  space  (or  veins)  as  a  function  of  pressure 
drop  across  the  air/liquid  interface  (i.e.,  capillary  ten¬ 
sion).  This  work  was  in  part  prompted  by  Mader’s  (1992 
a,  b)  theoretical  and  experimental  study  of  equilibrium 
vein  widths  in  temperature  ice  above  -0.5°C. 

EXPERIMENTAL  PROCEDURE 

We  measured  the  moisture  characteristic  curve  (MC) 
of  ice  with  a  procedure  similar  to  that  used  for  soils. 
Figure  1  shows  a  schematic  of  the  experimental  appa¬ 
ratus,  in  which  an  ice  sample  is  placed  within  a  her¬ 


metically  sealed  Tempe  cell  (Fig.  2)  and  then  subjected 
to  incremental  increases  in  air  pressure  to  “push”  off 
the  water.  To  control  the  ice  vein  width,  and  hence  the 
permeability  and  entrance  pressure,  we  submerged  the 
ice  sample  in  a  refrigerated  bath  maintained  just  below 
0°C.  The  setup  is  kept  within  a  coldroom  maintained 
near  +5°C.  The  Tempe  cell  consists  of  a  brass  cylinder 
and  two  plastic  end  pieces,  sealed  airtight  with  pairs  of 
“O”  rings.  Water  volume  removed  from  the  ice  is  mea¬ 
sured  through  level  changes  in  a  graduated  pipette  con¬ 
nected  by  plastic  tubing  to  the  bottom  of  the  Tempe 
cell.  A  1-bar  semipermeable  ceramic  plate  (Soilmoisture 
Equipment  Corp.)  separates  the  pressurized  air  cham¬ 
ber  from  the  open  water  reservoir.  Bubbling  pressure 
and  hydraulic  conductivity  of  the  plate  are  specified  as 
20  to  30  PSI  and  3.46  x  10“^  cm/sec.  Water  remained 
unfrozen  within  the  ceramic  stone  because  of  its  relatively 
high  matric  tension.  A  resistance-wire  heater,  connected 
to  a  temperature  regulator,  prevented  freezeup  of  water 
within  the  submerged  tubing. 

We  used  a  custom-made  refrigerated  bath  for  our 
experiments  (Black  1986),  consisting  of  a  20-gallon 
container  that  is  subjected  to  continual  cooling  by  a 
refrigeration  system.  Because  of  its  large  thermal  iner¬ 
tia,  this  bath  is  capable  of  constant  temperature  condi¬ 
tions  (±  0.003 °C)  for  time  durations  on  the  order  of 
weeks.  A  Bailey  temperature  controller  is  used  to  apply 
the  necessary  heat  through  an  immersed  heating  ele¬ 
ment  to  bring  the  bath  to  the  desired  temperature. 
Although  this  approach  is  not  energy  efficient,  it  works 
well  because  it  is  easier  to  precisely  warm  the  water 
than  it  is  to  cool  it.  To  promote  efficient  energy 
exchange,  the  bath  is  continually  mixed.  We  mounted 
the  Tempe  cell  apparatus  on  a  freestanding  support  to 
isolate  it  from  the  agitation  of  the  bath. 

We  made  experimental  runs  on  two  ice  samples.  The 
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Figure  2.  Ice  sample  is  placed  within  a  hermetically 
sealed  Tempe  cell. 


initial  run  tested  the  viability  of  the  experiment.  For 
this  run  we  used  a  quickly  prepared  ice  sample 
(described  in  the  next  section)  and  monitored  tempera¬ 
ture  of  the  ice  bath  with  a  single  thermistor  connected 
to  a  manual  readout.  High  sensitivity  of  water  levels  to 
bath  temperature  indicated  the  need  for  continual  and 
more  accurate  temperature  monitoring.  Thus,  on  the 
second  experiment  we  added  three  thermistors  and  a 
highly  accurate  (±  0.00 1°C)  platinum  resistance  ther¬ 
mometer  (PRT).  Data  from  the  thermistors  were  con¬ 
tinuously  collected  on  a  Campbell  CRIO  and  periodi¬ 
cally  transferred  to  a  computer.  The  PRT  was  manually 
monitored  with  a  S1220  Systemteknik  temperature 


meter.  The  full  combination  as  shown  in  Figure  1  was 
in  place  for  the  second  experiment. 

SAMPLE  PREPARATION,  INSTALLATION, 
AND  ANALYSIS 

We  prepared  two  ice  samples,  both  seeded  with 
frazzle  ice  nucleated  from  singly  distilled,  supercooled 
water.  To  grow  the  first  sample  (Sample  A),  we  placed 
a  beaker  of  frazil  slurry  in  a  brine  bath  maintained  at  a 
temperature  near  -20°C  in  a  coldroom.  Continuous 
agitation  with  a  magnetic  stirrer  removed  air  bubbles 
ahead  of  the  freezing  front.  This  procedure  produced  a 
heterogeneous  sample  of  coarse  to  fine-grained  ice  (Fig. 
3  a,  b).  Large  crystals  on  the  order  of  centimeters  grew 
in  from  the  sides  and  bottom  of  the  beaker  and  (at  a 
slower  rate)  at  the  top  of  the  slurry.  The  central,  fine¬ 
grained  portion  was  agitated  by  the  magnetic  stirrer, 
which  we  removed  for  the  final  freezing.  For  the  sec¬ 
ond  sample  we  sought  a  refined  technique  that  would 
produce  a  more  homogenous  and  fine-grained  sample. 
To  prepare  the  second  sample  (Sample  B),  we  placed  a 
beaker  of  frazil  ice  on  a  -25 °C  cold  plate  in  a  room 
where  the  temperature  was  kept  slightly  above  0°C. 
Water  froze  from  the  bottom  up  while  being  stirred  at  3 
rpm  with  a  motorized  paddle  held  just  above  the  freez¬ 
ing  front  on  an  adjustable  mount.  This  procedure  grew 
relatively  fine-grained  ice  (Fig.  3c-d)  at  a  rate  of  approx¬ 
imately  1  cm  per  hour.  Both  samples  were  disc-shaped. 
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a.  Bottom  of  Sample  A. 


b.  Top  of  Sample  A. 


c.  Bottom  of  Sample  B.  d.  Top  of  Sample  B. 

Figure  3.  Thin  sections  of  ice  samples  under 
polarized  light. 


about  8.3  cm  in  diameter  and  3.6  to  4.3  cm  high. 
Samples  A  and  B  weighed  180  and  222  grams,  respec¬ 
tively. 

To  assemble  the  Tempe  cell,  we  filled  the  pipette, 
tubing,  and  the  attached  bottom  plate  with  singly  dis¬ 
tilled  water  near  the  freezing  point.  We  then  moved  the 
apparatus  to  a  coldroom  to  add  the  ice  sample  and  secure 
the  remaining  parts  of  the  cell.  To  ensure  good  contact, 
we  melted  the  bottom  of  the  ice  sample  slightly  before 
placement  on  the  porous  stone.  During  installation, 
water  rose  to  fill  a  small  gap  between  the  ice  sample 
and  the  8.6-cm-diameter  brass  cylinder.  After  comple¬ 
tion  of  the  experiments,  ice  samples  were  returned  to 
the  coldroom  and  removed  from  the  Tempe  cell  for 
analysis. 

Figure  4  shows  photographs  of  the  rotted  ice  samples 
taken  immediately  after  removal  from  the  ice  bath.  Both 
samples  show  pronounced  grain-boundary  grooves  with 
widths  on  the  order  of  1  mm.  Grooves  initially  appeared 
more  pronounced  than  in  the  photographs,  perhaps 
because  refreezing  had  already  started  to  occur  before 
the  photographs  were  taken.  The  bubbly  central  area  in 
Sample  A  indicates  that  pressurized  air  entered  the  ice. 
Sample  B  was  somewhat  cloudy,  but  showed  no  dis¬ 
tinct  bubbly  regions.  Several  candled  grains  along  the 
sample  sides  flaked  off  with  handling  (Fig.  4c).  Sample 
B  also  evidenced  possible  surface  melting  towards  the 
top  of  one  side. 

We  prepared  thin  sections  to  analyze  the  ice  struc¬ 


ture.  Figures  3a-d  display  the  grain  cross  sections  under 
polarized  light.  Prior  to  sectioning,  the  ice  was  stored 
for  several  days  at  -20°C  and  thus  lost  most  evidence 
of  rotting.  By  counting  the  number  of  grain-boundary 
intercepts  per  length  of  randomly  oriented  lines  and 
taking  the  reciprocal,  we  determined  grain  sizes  (G)  of 
3.9  mm  and  3.2  mm  for  Samples  A  and  B,  respectively. 

We  melted  the  remains  of  the  ice  samples  and  meas¬ 
ured  the  impurities  concentration  with  a  Dionex  ion 
chromatograph.  This  analysis  yielded  a  chloride  con¬ 
centration  of  35.8  ppm  for  Sample  A  and  5.4  ppm  for 
Sample  B,  giving  respective  bulk  molalities  of  1.01  x 
10“3  moles  kg“l  and  1.52  x  10^  moles  kg“k  Both  con¬ 
centrations  were  higher  than  expected  for  singly  dis¬ 
tilled  water.  The  Hanover  Water  Company,  for  instance, 
reported  C\~  concentrations  of  around  0.35  ppm  for  tap 
water  at  the  fire  station  near  CRREL  and  Mader  ( 1 992b) 
measured  concentrations  two  to  three  orders  of  magni¬ 
tude  less  for  her  laboratory-grown  samples.  Mader, 


a.  Top-down  view  of  Sample  A. 


b.  Top-down  view  of  Sample 
B  within  the  Tempe  cell. 


c.  Side  view  of  Sample  B  with 
its  top  downward.  Note  candled 
ice  grains  and  possible  surface 
melting  at  lower  right. 

Figure  4.  Rotted  ice  samples  from 
dismantled  experiment. 
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Table  1.  Geometric  definitions  and  relationships  of  ice  structure.  All  units  are  SI. 

Parameter 

Functional  relationship 

Evaluated 

G  (grain  size) 

d)  (dihedral  angle) 

length/grain  boundary  intercepts 

32° 

0  (contact  angle) 

0° 

Gji ,  Oai ,  and  Giv  (surface  tension) 

0.034,  0.075,0.109  (Jm-2) 

VEINS 

a 

a  =  7i/6-^/2 

14° 

dy  (vein  width) 

Ty  (radius  of  curvature  of  vein) 

ry  =  dy/2  sina 

ry  =  (7y/0.4838  (m) 

Ay  (cross  sectional  area  of  vein) 

\  sin^  a  +  -^sin2a  -  3aj 

=  0.0725  rA  (m2) 

/  (length  of  veins  per  unit  volume)^ 

l  =  3/(fi 

Rii  (hydraulic  radius) 

2  vein  area/vein  circumference 

/C  - 

Ru  =  0.2044  Jv(m) 

k[  (permeability) 

‘  2-j2G^(\  +  b^) 

tav  (vein  air-entry  radius) 

r^v  =  (<^v  !  2)  71  /  6  -  tan  a  /  2) 

tav  =  0.2273  dy  (m) 

NODE 

(nodal  volume/ vein  volume)  ^ 

n.^9dylG 

rAn  (node  air-entry  radius)^ 

rA„=«/2)^3-8.61)/(4jt) 

rAn  =  0.6357  dy  (m) 

LENS 

Fl  (volume  of  lens)^ 

Vl  =  ■^|^2-3cos-^  +  cos^ 

0.0093  ri}  (m3) 

Sy  (surface  to  volume  ratio  of  ice  grain)^ 

(eVs +3)/(4V2g) 

2.37/G 

GROOVE  (Fig.  8c) 

P 

sin“l  [(rgCOs(d)/2)/(ra  +  rg)] 

Nye  (1991b),  ^Nye  (1989),  ^Nye  and  Mae  (1972),  dLliboutry  (1996). 

however,  discarded  the  dirtier  portions  of  her  samples 
before  analysis.  Contamination  from  the  preparatory 
brine  bath  used  to  grow  Sample  A  undoubtedly  con¬ 
tributed  to  the  high  concentration  of  C\~.  For  Sample 
A,  we  assume  that  the  C\~  cations  are  located  within 
the  ice  lattice  (Wolff  1996)  and  that  an  equal  number  of 
Na+  anions  are  apportioned  50%  each  between  grain 
boundaries  and  veins  (Harrison  and  Raymond  1976). 
Both  samples  also  accreted  impurities  during  post- 
experimental  handling  and  sectioning  (e.g.,  Fig.  4c). 
The  proportional  impact  of  this  contamination  is  greater 
for  Sample  B. 

Because  expelled  impurities  concentrate  ahead  of 
the  freezing  front,  they  occur  within  a  central  core  in 
Sample  A  and  towards  the  top  of  Sample  B.  Mader 
(1992b)  found  that  89%  of  the  impurities  resided  in 
14%  of  the  ice  volume.  With  these  proportions,  the  bulk 
molality  for  the  “clean”  part  of  the  sample  is  reduced 
to  a  factor  of  0.13  while  that  for  the  “dirty”  part  is 
enhanced  to  a  factor  of  6.4.  Here  we  assume  that  impu¬ 
rities  are  immobile,  although  diffusion  of  solutes  is  pos¬ 


sible  at  these  molality  levels.*  The  vein  bulk  molality 
Wy  for  the  clean  part  of  Sample  A  is  then  6.46  x  10“^ 
moles  kg“l .  Because  of  possible  post-experimental  con¬ 
tamination,  we  take  as  1.00  x  10“^  moles  kg“l  for 
the  clean  part  of  Sample  B,  in  line  with  Mader  (1992b) 
and  closer  to  the  value  computed  from  Hanover  tap 
water. 

THEORY 

Ice  structure  is  determined  thermodynamically  from 
combined  thermal  and  mechanical  stresses  in  the  ice. 
Interfacial  surface  tensions  define  the  geometry  of  air 
and  water  inclusions.  Table  1  summarizes  the  geometri- 


*  In  colder  ice,  the  diffusivity  of  heat  is  much  greater  than  the 
diffusivity  of  salt  and  limits  the  diffusion  of  impurities.  As  tempera¬ 
tures  approach  0°C  and  with  increasing  molality,  however,  latent  heat 
affects  in  ice  lower  its  thermal  diffusivity  (Nye  1991a).  Sample  A 
has  a  sufficiently  high  molality  that  thermal  and  salt  diffusivities  are 
of  the  same  magnitude  for  a  temperature  around  -0.0 1°C  and  thus 
some  diffusion  of  solutes  may  occur. 
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cal  definitions  and  relationships  used  in  this  study.  Liq¬ 
uid  inclusions  occur  as  veins  along  three-grain  bound¬ 
aries,  as  nodes  where  four  veins  join,  and  as  lenses  be¬ 
tween  two-grain  boundaries.  Veins  dilate  to  widths  of 
tenths  of  millimeters  as  ice  nears  0°C,  and  thus,  together 
with  nodes,  form  a  permeable  network  for  the  flow  of 
water.  Nye  (1989)  provides  a  detailed  description  of 
the  geometry  of  the  vein-node  system.  Veins  have  a 
concave  triangular  cross  section  (Fig.  5a)  which  is  fixed 
by  the  dihedral  angle  and  does  not  vary  with  vein  size. 
Nodes  take  the  shape  of  tetrahedrons  with  nonspherical 
faces  and  open  corners  (Fig.  5b).  Water  lenses  form 
spherical  caps  (Fig.  5c),  where  divergence  at  the  ice- 
water-ice  juncture  is  again  fixed  by  the  dihedral  angle. 
Because  of  the  respective  concave  and  convex  shapes 
of  the  ice/water  interfaces,  water  is  under  tension  in 
veins  and  under  pressure  in  lenses. 

The  phase-equilibrium  temperature  of  water  inclu¬ 
sions  is  defined  by  pressures  in  the  ice  and  liquid  phases 
and  by  impurities  in  the  water.  The  change  in  equilib¬ 
rium  temperature  AT  from  its  bulk  melt  value  Tq 


(=  273 . 1 5  kelvins  or  0°C)  computes  from  the  Clausius- 
Clapeyron  equation  as 


M  = 


T'o(p\-  Po 

k\  [  p\ 


Pj-Po^ 
Pi  y 


-Kf 


c 


(1) 


where  L[\  is  the  latent  heat  of  freezing  (3.335  x  lO^J 
kg“l),/>o  is  the  reference  pressure  (atmospheric), is 
the  water  pressure,  pi  is  the  ice  pressure,  p|  is  the  den¬ 
sity  of  water  (1000  kg  m~^),  pi  is  the  density  of  ice 
(917  kg  m“^),  C  is  the  molality  of  the  liquid  water  in¬ 
clusion  (moles  kg“l),  andiCp  is  the  modal  freezing  point 
depression  constant  for  water  (=  1.855  K  kg  mole“l). 
Thus,  increases  or  decreases  in  the  pressure  of  ice  or 
water,  respectively,  lower  the  equilibrium  temperature. 
The  Laplace  equation  defines  the  pressure  drop  across 
the  ice/water  interface  as pi-p\  =  2ai\  Ir  ii ,  where  r  ii  is 
the  mean  radius  of  curvature  of  the  interface  and  an  its 
surface  tension  (0.034  J  m~^).  Taking  the  excess  pres¬ 
sure  in  the  ice  as  the  gauge  pressure,  pQ,  and  defining 
p\  by  the  Laplace  equation,  eq  1  becomes: 


b.  Perspective  drawing  of 
node.  (After  Nye  1989.) 


Water 


c.  Lens  cross  section. 


Figure  5.  Geometry  of  water  inciusions  in  ice. 


AT  =  ^ 
Li\ 


PG 


vPl 


20i| 

Pi  'll 


-KpC. 


(2) 


The  terms  in  eq  2  represent,  respectively,  temperature 
depression  due  to  increased  gauge  or  “atmospheric” 
pressure,  interfacial  curvature,  and  impurities.  We  as¬ 
sume  that  the  gauge  pressure  acts  equilaterally  in  the 
ice,  but  note  that  this  is  not  necessarily  the  case. 

The  molality  Cy  of  the  vein-node  system  computes 
as  niy/vy,  where  niy  is  the  vein  bulk  molality  and  Vy  is 
the  fractional  volume  of  liquid.  Applying  thermody¬ 
namic  reasoning,  Nye  (1989)  computes  the  ratio  of  liq¬ 
uid  volume  in  nodes  to  that  in  veins  as  13.89(7y/G.  The 
fractional  volume  of  liquid  in  the  vein-node  system  is 
then^y/  (1  +  6.719ry/G),  whereby  is  the  cross-sectional 
area  of  a  vein  and  /  the  vein  length  per  unit  ice  volume. 
Using  the  geometric  relationships  in  Table  1,  eq  2,  and 
the  Laplace  equation,  the  equilibrium  temperature  for 
veins  is  (Lliboutry  1996,  Mader  1992b*) 


Ar^  =-7.413  X 10“^  j3g- 


2.785  xlQ-S 
'V 


8.528  my 
r2(G  + 6.7 19rJ 


(3) 


where  is  the  radius  of  curvature  of  the  vein.  Figure  6 


*  In  deriving  eq  3,  we  do  not  include  pumping  of  solutes  by  freeze- 
thaw  contraction  and  expansion  of  the  ice  volume.  This  and  smaller 
grain  sizes  for  our  samples  yield  smaller  vein  widths  than  observed 
by  Mader  (1992b). 
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Temperature  Depression  {°C) 

Figure  6.  Vein  diameter  vs.  equiiibrium  temperature  for  a  range  of  vein 
buik  moiaiities  rriy.  Computed  from  eq  3,  where  G  =  0.0032  m  and  po  =  0. 


shows  vein  diameter  dy  (=  0.4838  r^)  as  a  function  of 
temperature  for  a  range  of  molalities  that  bracket  our 
samples  and  Figure  7  shows  the  corresponding  volume 
fraction  of  liquid  water.  The  impurities  term  dominates 
the  curvature  term  below  about  -1  x  10“^  °C.  For  these 
plots,  the  gauge  pressure  (or  excess  air  pressure)  is  taken 
as  0.  At  the  maximum  gauge  pressure  of  about  1  bar  or 
15  PSI,  the  equilibrium  temperature  would  be  lowered 
overall  by  0.007°C  from  that  in  Figure  7.  We  note  that 
the  liquid  volume  fraction  is  quite  low  for  the  lower 
molalities  and  suggest  that  water  released  in  the  MC 
curves  was  primarily  supplied  by  the  contaminated 
portions  of  the  samples. 

The  molality  Cl  of  a  lens  is  7z(ri^sin<^/2)^mi^/ 


(0.5SyV]J,  where  Sy  is  the  surface-to-volume  ratio  of 
an  ice  grain,  tl  is  the  radius  of  the  lens  interface,  wl  is 
the  lens  bulk  molality ,  and  vl  is  the  lens  volume.  Again 
using  the  relationships  in  Table  1 ,  the  depression  tem¬ 
perature  for  a  lens  is 


AL  =-7.413x10-^0- 


5.570x10-^ 

>L 


40.19GmL 


(4) 


Although  the  temperature  differences  between  veins  and 
lenses  are  small,  the  resulting  heat  fluxes  have  impor¬ 
tant  consequences  for  the  ice  structure.  If  the  gauge 


Temperature  Depression  (°C) 

Figure  7.  Volume  fraction  of  liquid  water  vs.  equilibrium  temperature  for  a 
range  of  vein  bulk  molalities  mv 
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pressure  acts  equilaterally  within  the  ice  and  the  impu¬ 
rities  concentrations  are  similar,  then  heat  is  directed 
towards  the  veins  and  the  lenses  will  disappear.  Impu¬ 
rities  in  the  veins  are  potentially  mobile,  however,  and 
thus  wl  may  exceed  in  the  contaminated  region  of 
the  sample.  Heat  is  then  directed  toward  the  colder 
lenses  and  the  veins  will  shrink.  If  wl  =  2my  and  = 
100  pm,  the  equilibrium  lens  radius  is  13  pm.  Thus 
while  lenses  may  form,  it  is  unlikely  that  they  will  grow 
large  enough  to  coalesce.  We  must  also  consider  the 
possibility  that  the  stress  does  not  act  equilaterally.  Nye 
and  Mae  (1972)  observed  the  formation  of  lenses  and 
shrinkage  of  veins  following  the  short  application  of  a 
compressive  stress  of  1.3  bar.  They  theoretically  com¬ 
pute  a  critical  overburden  pressure  of  0.23  bars  for  this 
effect  to  occur. 

In  this  paper  we  investigate  unsaturated  ice  and,  thus, 
must  also  take  into  the  account  the  thermodynamic  state 
of  ice  when  the  vapor  phase  is  present.  In  many  ways, 
the  saturated  and  unsaturated  regimes  in  ice  are  analo¬ 
gous  to  the  funicular  and  pendular  regimes  in  snow 
(Colbeck  1973,  Morris  1993).  Here  we  conceptualize 
ice  as  a  porous  medium,  similar  to  snow  or  soil,  where 
the  pores  are  either  vein  segments  or  nodes.  Each  pore 
has  a  critical  entry  pressure  that  must  be  reached  before 
air  can  replace  water.  For  a  cylindrical  pore  of  radius 
R,  the  radius  of  curvature  of  the  air/liquid  interface  is 
7?/cos0,  where  0  is  the  contact  angle  (Fig.  8a).  Without 
defining  in  detail  the  geometry  of  the  interfacial  sur¬ 
face  in  a  vein,  we  specify  an  equivalent  hemisphere 
whose  radius  r^v  is  that  of  an  inscribed  circle  (Fig.  8b). 
Although  veins  have  a  flared,  trumpet-like  shape  at  the 
surface  (Nye  1991b),  here  we  consider  only  interior 
veins  with  parallel  faces.  If  the  contact  angle  is  0,  the 
pressure  drop  across  the  interface  is pg,-p\  =  2aai  /^av, 
where  an  is  the  surface  tension  (0.075  J  m~^)  and  tav  = 
0.2273  dy  (Table  1).  This  pressure  drop  defines  the  air- 
entry  pressure  T^av 

PAY  =  0M00/dy  (5a) 

where  Pj^y  is  in  Pascals  and  dy  in  meters.  A  column  of 
air  replaces  water  in  the  vein  when  the  pressure  drop 
exceeds  Pj^y.  We  note  that  air  columns  of  smaller  radius 
not  touching  the  ice  faces  are  mechanically  unstable. 
The  newly  formed  cylindrical  meniscus  is  of  higher 
tension  than  the  hemisphere  and  immediately  retreats 
into  the  three  grain-boundary  vertices,  as  shown  in  Fig¬ 
ure  8b.  Water  inclusions  in  the  unsaturated  regime  thus 
consist  of  concave  triangular  columns  in  the  grain¬ 
boundary  grooves.  Water  cannot  refill  the  vein  until  the 
gauge  pressure  is  lowered  to  the  point  Pay^2  where  the 
menisci  again  touch.  Thus  the  water  re-entry  pressure 
is  half  the  air-entry  value: 


R  (cosO)"^ 


R 


a.  Air/liquid  interface  within  a  cyiinder. 


b.  Air-entry  radius  of  curvature  with¬ 
in  a  vein  and  hypothesized  iocation 
of  iiquid  in  an  unsaturated  vein. 


c.  Grain-boundary  groove. 


Figure  8.  Geometry  of  iiquid  inciusions  in  contact 
with  air  (not  to  scaie). 

P^y  =  0.3300/ dy.  (5b) 

The  configuration  shown  in  Figure  8b  is  not  in  mechani¬ 
cal  equilibrium.  We  discuss  potential  changes  in  ice 
curvature  in  the  following  paragraph.  Similarly,  we  use 
the  radius  of  the  inscribed  sphere  within  a  node  to  com¬ 
pute  its  air-entry  pressure  T^An-  Using  the  nodal  geom¬ 
etry  determined  by  Nye  (1989), 
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0.2359/ Jv  (5c) 

roughly  one  third  that  of  Pj^y.  The  water  re-entry  pres¬ 
sure  is  again  half  of  this. 

P^^=0.m0/dy.  (5d) 

In  a  series  of  interconnected  pores,  the  smallest  pore 
will  block  the  entry  of  air  and  the  largest  pore  will  block 
the  re-entry  of  water.  In  this  idealized  scheme  of  alter¬ 
nating  veins  and  nodes,  the  air-entry  pressure  and  water 
re-entry  would  be  P\y  =  0.6600/dy  and  T*wn  =  0. 1 1 80/ 
dy.  For  a  vein  diameter  of  0. 1  mm  these  evaluate  as  6.6 
kPa  and  1.2  kPa.  This  scheme  thus  predicts  a  highly 
hysteretic  moisture  tension  curve  for  ice. 

When  air  replaces  water  within  the  vein  system, 
mechanical  and  thermal  equilibrium  at  the  newly 
formed  triple-phase  junction  is  suddenly  altered  (point 
A  in  Fig.  8b).  Application  of  the  Laplace  equations  to 
the  three  interfaces  defines  the  mechanical  equilibrium 
(Colbeck  1979) 


where  Gia  is  the  ice-air  surface  tension  (0.109  J  m~^), 
and  r  ai ,  r  {\ ,  and  ria  are  the  respective  mean  radii  of 
curvature  of  the  three  interfaces.  Figure  8c  shows  an 
open  groove  configuration,  where  pressurized  air 
bounds  the  upper  liquid  surface.  Here  rg\  =  2r^,  ^  i\  = 
2rg  and  =  2ry  ,  where  is  the  radius  of  curvature 
of  the  exposed  ice  surface.  The  higher  surface  tension 
of  the  air/liquid  interface  dictates  the  groove’s  flared 
shape.  Within  a  vein,  these  grooves  must  configure 
within  an  equilateral  triangle,  as  shown  in  Figure  8b. 
In  order  to  meet  this  additional  restraint,  we  derive  a 
somewhat  surprising  solution  (not  included  in  the  re¬ 
port)  in  which  the  curvature  of  the  vein  walls  changes 
from  positive  to  negative.*  If  pressures  in  the  liquid 
and  ice  are  expressed  relative  to  gauge  pressure,  eq  1 
defines  a  new  equilibrium  temperature  for  the  ice/wa¬ 
ter  interface: 


where  Cg  is  the  concentration  of  impurities  in  the  groove 
water.  For  a  similar  impurities  concentration,  the  tighter 
curvature  and  increased  surface  tension  of  the  air/liq¬ 
uid  interface  imply  a  lower  equilibrium  temperature 
relative  to  the  saturated  state  (i.e.,  eq  3).  Heat  flow 
directed  towards  this  colder  region  fosters  melt. 
Mechanical  equilibrium  and  increased  pressure  also 
suggest  surface  melting  of  the  newly  exposed  vein  wall. 
Without  proposing  a  detailed  equilibrium  geometry  for 
the  unsaturated  vein,  we  raise  the  possibility  that  veins 
exposed  to  air  will  change  curvature  and  enlarge.  This 
process  is  limited  by  dilution  of  the  impurities  by  melt¬ 
water. 

To  describe  water  flow  in  ice,  we  also  need  to  esti¬ 
mate  the  permeability  of  the  vein  system.  We  use  the 
hydraulic  radius  rn  (2  vein  area/vein  circumference)  to 
describe  an  equivalent  radius  for  a  cylindrical  conduit. 
In  terms  of  vein  geometry  (Fig.  5a),  th  =  0.2044  dy. 
Lliboutry  (1996)  computes  permeability  from 
Pouseilles’  equation,  multiplied  by  the  vein  density  and 
divided  by  a  tortuosity  factor.  Since  water  flow  is  con¬ 
stricted  by  the  smallest  pores  in  the  network,  we  also 
add  a  vein  diameter  distribution  factor  (Dullien  1992). 
If  the  vein  system  consists  of  vein-node  doublets  with 
sizes  dy  and  dylb,  the  harmonic  mean  permeability  is  2/ 
(1  +  b^)  times  that  with  the  larger  size  alone.  Taking 
vein  density  as  2IG'^  and  tortuosity  as  V2  ,  the  ice  per¬ 
meability  is  then 

_  1.94x10-3^/4 

‘  ^  2a/2G2(1  +  />4)“  G^{\  +  b*) 

Thus  for  dy  =  5  X  10“^  m,  G  =  3.5  x  10“^  m,  and  b  =  2, 
is  5.82  X  10“!^  rrP.  The  permeability  of  ice  is  below 
that  for  soils  of  comparable  pore  size  because  of  the 
low  vein  density. 

ki  computed  with  eq  9  is  the  intrinsic  or  water-saturated 
value.  Air  blockages  in  unsaturated  ice  can  greatly 
reduce  this  value. 

EXPERIMENTAL  RUNS 


AT  = 


- 


Ai 


r  1 

n 

^al 

^ia 

Pg 

^pi 

Pi  J 

Pi '■a 

PiC_ 

-K^q  (7) 


which  evaluates  as 


A^g  = -7.413x10-^0  - 


6.143x10-8 


9.736x10-8 


- 1.855  Cg 


(8) 


*  The  solution  is  available  from  the  authors  upon  request. 


We  made  several  runs  for  each  of  the  samples,  using 
the  experimental  setup  described  above.  Several  of  these 
were  marred  by  temperature  fluctuations  and  are  not 
reported  here.  For  each  run,  we  incrementally  cycled 
the  gauge  pressure  from  0  through  1  bar  and  back.  Each 
pressure  level  was  held  for  about  two  days  or  until  the 
water  level  stabilized.  We  report  only  the  equilibrium 
water  level.  Because  of  the  long  equilibrium  times  for 
ice,  each  run  required  many  weeks  to  complete.  Prior 
to  a  run,  we  rotted  the  ice  sample  for  several  days  at  a 
fixed  bath  temperature  to  obtain  a  given  vein  dilation, 
as  shown  in  Figure  6.  We  noted  the  amount  of  water 
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0  12  3  4 

Water  Removed  (cm^) 

Figure  9.  Volume  of  water  exiting  (open  triangles)  and  re-entering 
(solid  triangles)  Sample  A  vs.  gauge  pressure  (pq)  for  bath 
temperatures  between  -0.09°C  and  -0.08°C.  The  bath  temperature 
drifted  upwards  0.01  to  0.02°C  at  pressure  levels  of  20  and  47  kPa  on 
the  drying  curve. 


forced  from  the  ice  by  changes  of  water  level  in  the 
pipette.  Temperature,  pressure  level,  and  water  level 
were  monitored  one  to  three  times  a  day  and  tempera¬ 
ture  continuously  for  Sample  B. 

Figure  9  shows  results  from  the  initial  run  on  Sample 
A  with  the  bath  temperature  between  -0.09  and  -0.08°C. 
The  plot  shows  that  ice  desaturates  and  rewets  in 
response  to  changes  in  incremental  pressure  in  a  fash¬ 
ion  similar  to  unsaturated  soil.  Temperature  fluctuations 
of  ±  0.0 1  °C  and  changes  in  gauge  pressure  caused  struc¬ 
tural  changes  in  the  ice  over  the  run.  In  addition,  the 
sample  was  not  homogenous.  Localized  impurities, 
grain  size  variations,  and  temperature  gradients  in  the 
bath  caused  a  potentially  wide  range  of  vein  sizes. 
Measured  temperature  gradients  in  the  bath  suggest  that 
the  top  of  the  sample  was  0.01  to  0.02°C  warmer  than 
the  bottom.  Because  of  rapid  conduction  through  the 
brass  cylinder,  sides  of  the  sample  were  the  most  sensi¬ 
tive  to  temperature  fluctuations.  We  hypothesize  that 
smaller  vein  sizes  towards  the  bottom  of  the  sample 
control  the  entrance  pressure  while  larger  veins  and 
nodes  in  warmer  and  contaminated  regions  control  the 
liquid  volume. 

Figure  9  shows  that  air  first  entered  Sample  A  at  a 
pressure  of  5  kPa  and  water  continued  to  be  pulled  from 
the  ice  up  to  a  pressure  of  57  kPa.  Again,  we  theoreti¬ 
cally  relate  the  vein  air-entrance  pressure  to  its  width 
as  Pj^y  =  0.6600/dy,  which  suggests  a  range  of  vein  sizes 
of  132  to  12  jam.  Consistent  with  this,  the  Clausius- 
Clapeyron  equation  (eq  3)  predicts  a  vein  width  of  127 
jam  for  the  clean  part  of  Sample  A.  Our  tenfold  range 
in  vein  widths  is  much  larger  than  the  twofold  varia¬ 
tion  observed  by  Mader  (1992a),  but  possibly  reflects 


sample  heterogeneity  and  temperature  gradients  in  the 
bath.  A  total  of  2.9  cm^  of  liquid  was  removed  from  the 
sample.  The  volume  fraction  Vy  for  clean  ice  (6.46  x 
10“^  moles  kg“l)  is  around  0.001  and  thus  cannot 
account  for  this  liquid.  A  volume  fraction  of  0.07  for 
the  contaminated  portion  (3.23  x  10“^  moles  kg“l) 
yields  2.2  cm^  of  water.*  Larger  increments  in  water 
at  gauge  pressures  of  22  and  49  kPa  correspond  to 
periods  of  upward  temperature  drift  of  0.01  to  0.02°C. 
After  these  occurrences,  we  reset  the  bath  temperature 
to  -0.09°C.  In  the  reverse  direction.  Figure  9  shows 
water  first  re-entering  the  sample  at  13  kPa.  Recall  that 
our  theoretical  estimate  of  water  re-entry  pressure  is 
T^Wn  ^0.11 80/dy,  which  predicts  that  water  will  re-enter 
the  12-|Lim  veins  at  a  pressure  of  10  kPa.  Figure  9  thus 
supports  our  hypothesis  that  the  MC  curve  of  ice  is 
highly  hysteretic.  Re-entry  of  the  remaining  1 .6  cm^  of 
water  at  a  pressure  of  5  kPa  occurred  slowly  over  a 
two-week  period,  suggesting  a  low  permeability  for  the 
ice  and  possibly  a  hysteresis  effect.  With  gauge  pres¬ 
sure  reduced  to  0,  an  additional  0.5  cm^  of  water  entered 
the  sample  under  the  2  kPa  gravitational  head  over  the 
next  two  days. 

During  the  next  two  weeks  we  again  increased  the 
gauge  pressure  up  to  1  bar.  Surprisingly,  the  water  flow 
direction  was  unaffected  by  this  reverse  in  pressure  and 
continued  to  drop  slowly  at  around  0.05  cm^  per  day. 
We  examine  this  anomalous  behavior  and  issues  con¬ 
cerning  ice  permeability  in  the  next  section. 

The  results  of  this  first  experiment  were  promising. 


*  Computed  as  0.07  x  222  cm^  (total  volume)  x  0. 14  (fraction  of  dirty 
ice). 
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Figure  10.  Volume  of  water  exiting  (open  triangles)  and  re-entering  (solid 
triangles)  Sample  B  vs.  gauge  pressure  (pg)  for  bath  temperatures  between 
-0.03°C  and  -0.01  °C.  The  bath  temperature  warmed  to  0.005°C  for  approximately 
two  hours  at  a  pressure  level  of  61  kPa. 

We  thus  endeavored  to  refine  the  experiment  by  creat¬ 
ing  a  purer  and  more  homogenous  sample  and  by  moni¬ 
toring  the  bath  temperature  more  closely.  For  the  sec¬ 
ond  experiment,  we  rotted  Sample  B  for  one  month  at 
-0.05°C  to  -0.06°C  and  then  began  the  first  of  four 
pressure  cycles.  None  of  these  were  particularly  suc¬ 
cessful.  In  the  first  run,  we  incremented  gauge  pres¬ 
sure  to  83  kPa  and  back  with  only  minor  changes  in 
water  level.  On  the  drying  curve  the  water  level 
increased  by  0.05  cm^  at  a  pressure  of  34  kPa  and  on 
return  curve  it  decreased  by  about  half  this  amount  when 
the  pressure  was  reduced.  The  corresponding  vein  width 
of  20  |Lim  is  consistent  with  a  molality  of  1  X  10“^  moles 
kg“l .  The  tolerances  in  water  level  are  so  small,  how¬ 
ever,  that  we  cannot  attribute  much  significance  to  this 
value. 

For  the  second  run  we  increased  the  bath  tempera¬ 
ture  to  -0.03°C  to  obtain  greater  vein  dilation.  Figure 
10  shows  the  results  of  this  run.  For  the  first  pressure 
increments  the  water  level  decreased  slightly,  reminis¬ 
cent  of  this  anomalous  occurrence  for  the  second  run 
on  Sample  A.  The  4.3-cc  rise  in  water  level  at  a  pres¬ 
sure  of  60  kPa  occurred  when  the  bath  temperature 
approached  0°C  and  slightly  exceeded  it  for  about  an 
hour.  At  the  time  of  observation  and  during  most  of  the 
elevated  period,  the  temperature  was  around  -0.0 1°C. 

Because  the  equilibrium  temperature  for  4.3  cm^  of  vein 
water  is  around  -0.00 1  °C  for  the  dirtier  part  of  sample, 
surface  melt  undoubtedly  contributed  to  this  volume. 

After  returning  the  bath  temperature  to  -0.03°C,  we 
incremented  the  pressure  to  94  kPa  with  no  further  loss 


of  water.  On  the  return  phase,  the  water  level  decreased 
overall  by  0.2  cm^  as  the  gauge  pressure  was  decreased 
to  4  kPa.  The  first  drop  at  80  kPa  corresponds  to  a  very 
small  grain  diameter  of  2  jam.  If  the  spread  in  vein  di¬ 
ameters  is  again  tenfold,  this  implies  a  maximum  size 
of  around  20  jam,  consistent  with  run  1 .  With  the  gauge 
pressure  at  0,  water  slowly  re-entered  the  sample  over 
the  next  month. 

During  the  last  week  of  the  previous  run,  the  bath 
temperature  had  warmed  to  near  -0.0 1°C.  Thus  when 
we  increased  gauge  pressure  to  10  kPa  to  begin  run 
three,  we  immediately  removed  over  9  cm^  of  meltwa¬ 
ter  in  just  a  few  hours.  We  reduced  the  pressure  to  0 
and  all  but  2.8  cm^  again  re-entered  the  sample  over 
the  next  month. 

For  the  final  run,  the  bath  temperature  was  mostly 
between  -0.02°C  and  -0.0 1°C.  About  2  cm^  of  water 
was  removed  at  a  gauge  pressure  of  17  kPa.  An  addi¬ 
tional  7  cm^  were  removed  as  we  incremented  the  gauge 
pressure  to  93  kPa.  Because  of  the  sizable  amounts  of 
water,  surface  melting  most  likely  contributed  to  the 
water  volume.  This  close  to  the  0°C,  pressure  melting 
becomes  significant  and  could  account  for  additional 
meltwater  with  increasing  gauge  pressure.  No  water  re¬ 
entered  the  sample  on  the  depressurizing  cycle  or  dur¬ 
ing  the  following  week  at  0  gauge  pressure. 

PERMEABILITY 

Data  points  in  Figure  1 1  show  a  time  series  of  water 
re-entry  volumes  for  the  final  months  of  runs  2  and  3 
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Figure  11.  Days  vs.  volume  of  water  re-entering  Sample  B.  Solid  circles  are  for 
run  2  and  open  circles  for  run  3.  Lines  show  computed  flow  rates  for  the  porous 
plate  alone  and  for  hypothetical  configurations  of  the  ice  sample. 


on  Sample  B.  Gauge  pressure  had  been  reduced  to  zero, 
but  the  height  of  water  in  the  pipette  created  a  pressure 
head  at  the  bottom  of  the  sample.  In  both  data  plots, 
water  re-entry  is  at  first  quite  slow,  rises  quickly  for  the 
next  week  or  so,  and  then  returns  to  a  slower  rate. 
Although  it  is  possible  that  melt-induced  structural 
changes  caused  this  flow  variation,  bath  temperatures 
were  fairly  steady  over  both  time  periods.  Because  the 
data  plots  are  smooth  and  consistent,  we  attempt  to 
model  the  water  flow  and  thereby  gain  insight  into  the 
ice  permeability.  Figure  1 1  shows  computed  flow  rates 
for  the  porous  plate  alone  and  for  hypothetical  configu¬ 
rations  of  the  ice  sample. 

Figure  12  shows  the  simplified  structure  we  assume 
for  the  analysis:  a  cylindrical  stackup  of  the  porous  plate 
and  a  slightly  narrower  ice  core,  subject  to  a  pressure 
head  at  the  bottom  of  the  porous  plate.  The  initial  gravi¬ 
tational  head  is  a  height  Hq  of  water  in  the  pipette.  The 
receding  water  level  z'  in  the  pipette  relates  to  the  ris¬ 
ing  level  z  in  the  sample  as  z'  =  (^Aj.  z,  where  Aj.  is  the 
ratio  of  the  sample  area  to  the  pipette  area  and  (|)  is  the 
ice  porosity  ahead  of  the  water  front.  We  assume  a  one¬ 
dimensional  treatment,  but  note  that  this  is  valid  only 
in  the  absence  of  a  radial  pressure  gradient  beneath  the 
plate.  A  uniform  pressure  head  is  possible  if  the  resis¬ 
tance-wire  heater  (Fig.  1)  kept  the  area  beneath  the  plate 
from  totally  freezing.  If  this  area  were  instead  blocked 
by  ice,  we  estimate  that  the  radial  gradient  would  slow 
the  upward  velocity  of  the  water  front  by  roughly  one- 
seventh  (e.g.,  Muskat  1937). 

We  model  the  water  flow  with  Darcy’s  law: 


Figure  12.  Idealized  ice  configuration  used  in  flow 
computations.  Not  to  scale. 
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(|)Az 

dt 


{z  +  Zf)kik^ 

k^z  +  hz^ 


(Vp  +  pig) 


(10) 


where  |LL  is  the  viscosity  of  water  (1.789  X  10“^),  g  the 
gravitational  constant  (9.8  m  S“2),  p|  the  density  of  water 
(1000  kg  m“3),  and  Vp  the  pressure  gradient.  The  quan¬ 
tity  in  the  brackets  is  the  mean  harmonic  permeability, 
where  is  the  ice  permeability,  Zp  the  thickness  of  the 
porous  plate,  and  the  plate  permeability.  Pressures 
at  the  top  of  the  rise  z  and  at  the  bottom  of  the  porous 
plate  are  respectively  the  nodal  entry  tension,  -  T^wn? 
and  the  gravitational  pressure,  p|  g(/zo  -  z').  These  act 
over  the  combined  distance  +  z.  Integrating  eq  10 
we  obtain  an  implicit  equation  of  rise  height  in  terms 
of  time 


z 


B  . 


= t 


(11) 


where 

P<t) 


Po+?\gckit^lk^  ^j2) 

and  where  ^  ^Wn  +  PI  g^O  ^  c  =  (I  +  (|)Aj.).  With 
varying  parameters  and  boundary  conditions  (Table  2), 
eq  11  describes  the  three  flow  phases  observed  in  Fig¬ 
ure  1 1 .  We  explain  the  three  phases  as  follows:  1)  water 
rises  under  tension  within  the  central  ice  core  until  the 
veins  are  saturated,  2)  water  rises  under  a  gravitational 
head  to  fill  the  air  gap,  and  3)  water  flows  to  a  seepage 
face  at  the  top  of  the  ice  sample. 

The  hypothesis  of  unsaturated  flow  in  phase  1  is  the 
most  speculative  of  the  three  scenarios.  Flow  under  ten¬ 
sion  prevents  entry  of  water  into  the  side  air  gap  and 


hypothetically  explains  the  observed  delays  of  approxi¬ 
mately  three  and  seven  days  in  Figure  11.  We  recall 
that  slow  equilibrium  times  also  prevented  complete 
resaturation  on  the  wetting  cycle  for  runs  1  of  both 
samples.  In  these  cases,  water  continued  to  enter  the 
ice  after  the  next  drying  cycle  was  initiated  (e.g.,  dry¬ 
ing  curve  in  Fig.  10).  Solution  of  eq  11,  subject  to  the 
observed  delay  times  and  a  rise  height  z  of  h[,  yields 
extremely  low  ice  permeabilities  of  around  4  X  10“l^ 
and  1  X  10“!^,  respectively,  for  runs  2  and  3.  Corre¬ 
sponding  air  entry  values  are  4.3  and  3.4  kPa.  Appen¬ 
dix  A  provides  details  on  the  solution. 

Once  the  sample  is  saturated,  water  rises  rapidly  to 
fill  the  air  gap.  Tension  at  the  top  of  the  rise  (-Pwn)  is 
now  0  and  the  porosity  is  1 .  We  determine  the  thick¬ 
ness  of  the  air  gap  tg  from  the  observed  volumes  of 
water  re-entry  during  phase  2.  Note  that  the  area  ratio 
Aj;  is  now  close  to  1,  since  we  limit  flow  to  the  air  gap. 
Solution  of  eq  11  for  a  rise  height  of  h[  and  for  the 
observed  times  to  complete  phase  2  yields 
permeabilities  of  1  X  10“l^  and  3  X  10“l^  m^.  Thus,  the 
predicted  permeability  of  the  newly  frozen  ice  in  the 
air  gap  is  much  larger  than  that  for  the  central  core. 

In  order  to  not  “over-fit”  the  data,  we  use  the  same 
permeabilities  for  phase  3  that  were  determined  for  the 
central  core  in  phase  1.  These  ki  values  slightly 
overpredict  the  flow  rates  for  runs  2  and  3.  With  b  =  4, 
a  middle  ki  value  of  7  X  10“l^  rrP  corresponds  to  vein 
doublet  widths  of  3 1  and  8  qm.  These  widths  are  in 
line  with  values  suggested  by  the  MC  runs. 

The  higher  permeabilities  for  phase  2  are  consistent 
with  greater  vein  dilation  in  ice  at  the  sides  of  the 
sample,  but  also  suggest  that  by  comparison  ice  in  the 
central  core  underwent  structural  changes  during  the 
experiments.  All  four  runs  with  sufficient  data  for  analy¬ 
sis  showed  long  equilibrium  times,  on  the  order  of 
weeks,  for  the  wetting  curve.  We  speculate  that  vein- 
wall  melting  during  the  desaturation  process  may  have 


Table  2.  Parameters  used  in  computing  flow  rates  through  Sample  B. 


Phase  1 

Phase  2 

Phase  3 

hp  (m) 

0.00635 

h\  (m) 

0.043 

Rc(m) 

0.043 

kp  (m2) 

6.35X  lO-^e 

ho  (m) 

0.25  0.29 

0.25  0.29 

0.21  0.25 

Ar 

58.1 

0.837  0.884 

58.1 

tg  (m) 

0.0003 

0.0003 

0 

k\ 

4x10-''9  IxIO-^s 

1x10-''5  SxIO-^s 

4x10-''9  IxIO-'is 

b 

4.0  4.0 

na 

na 

<I> 

0.001  0.001 

1.0 

1.0 

P wn  (kPs) 

4.3  3.4 

0 

0 
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led  to  shrinkage  of  nearby  veins,  thus  increasing  the 
spread  of  vein  width  sizes.  Blockage  of  water  flow  by 
entrapped  air  is  another  likely  cause  of  reduced  ice  per¬ 
meability.  Both  of  the  causes  presume  that  air  enters 
the  veins.  Redistribution  of  impurities  would  lead  to 
vein  shrinkage  and  growth,  respectively,  in  contami¬ 
nated  and  purer  regions  of  the  samples. 

CONCLUSIONS 

From  the  relative  success  of  the  MC  experiment  and 
from  physical  evidence  of  a  large,  bubbly  region  within 
the  rotted  ice,  we  conclude  that  air  did  permeate  Sample 
A.  Although  water  volumes  were  somewhat  affected 
by  temperature  fluctuations,  our  data  suggest  that  ice 
desaturates  and  rewets  under  pressure  in  a  fashion  simi¬ 
lar  to  unsaturated  soil.  Air  entry  and  water  re-entry  pres¬ 
sures  generally  supported  our  theoretical  estimates  (eq 
5a  and  5d)  and  predictions  from  the  Clausius-Clapeyron 
equation  (eq  6).  The  observed  MC  curve  for  this  one 
experiment  was  highly  hysteretic,  in  line  with  our  ide¬ 
alized  flow  path  model  of  alternating  vein  segments 
and  nodes.  Because  of  small  water  volumes  and  inter¬ 
vening  melt  episodes,  MC  results  were  much  less  con¬ 
clusive  for  Sample  B.  The  few  usable  data  points,  how¬ 
ever,  are  again  in  line  with  theoretical  estimates.  Ancil¬ 
lary  information  on  flow  rates  through  Sample  B  cor¬ 
roborates  the  MC  runs  and  supports  the  hypothesis  that 
air  entered  the  sample.  For  two  runs,  we  speculate  that 
water  re-entry  into  the  side  air  gap  was  delayed  by  three 
and  seven  days  until  the  sample  was  resaturated.  For 
another,  water  levels  continued  to  drop  when  pressure 
was  reversed  upwards,  thus  suggesting  the  continued 
flow  of  water  under  tension  into  unfilled  airways. 

We  propose  that  water  inclusions  reside  in  the  verti¬ 
ces  of  air-filled  veins  and  take  the  form  of  concave  tri¬ 
angular  columns.  Thermal  and  mechanical  equilibrium 
requires  a  lowered  equilibrium  temperature  and  perhaps 
a  reversal  in  curvature  of  the  vein  wall  through  melt¬ 
ing.  The  formation  of  these  near-cylindrical  conduits 
has  important  implications  for  the  permeability  of  ice 
in  non-laboratory  settings.  Finer-grained  soils  under¬ 
lying  basal  ice  could  draw  water  from  the  veins  and,  in 
warm  spring  conditions,  began  this  process  of  vein 
enlargement.  Our  laboratory  results  show  that 
resaturated  ice  under  hydrostatic  tension  is  nearly  imper¬ 
meable,  with  permeabilities  on  the  order  of  1  X  10“l^ 
nfi.  Low  permeabilities  may  result  from  air  entrapment. 

We  conclude  that  moisture-tension  experiments  on 
ice  are  doable  and  offer  an  alternative  tool  for  studying 
ice  thermodynamics  in  the  temperate  regime.  Tempera¬ 
ture  fluctuations  were  the  major  drawback  in  this  experi¬ 
ment  and  could  be  stabilized  by  adding  an  encapsulat¬ 
ing  brine  bath  around  the  Tempe  cell.  We  suggest  in¬ 
creasing  the  diameter  of  the  sample  and  Tempe  cell  to 


gain  more  water  outflow  and  to  use  a  pipette  with  finer 
gradations.  Controlled  doping  of  the  sample,  unwit¬ 
tingly  accomplished  for  Sample  A,  would  also  increase 
the  water  volume.  Growing  a  larger  sample  and  dis¬ 
carding  the  contaminated  ice  would  ensure  greater  ho¬ 
mogeneity. 
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APPENDIX  A:  SOLUTION  OF  EQUATION  12  FOR  FLOW  UNDER  TENSION 


We  model  the  water  pathway  with  the  idealized  vein  doublet  discussed  above,  where  the 
dimensions  of  the  smaller  vein-node  pair  are  l/Z^th  the  larger.  Repeating  eq  9: 


_  1,94x10-3^/4 
G2(1  +  Z?4) 


(Al) 


where  G  =  3.2  X  10  The  volume  v  of  water  re-entering  the  sample  is  where  V 
(=  250  cc)  is  the  sample  volume.  With  v  measured  as  0.2  cc  we  take  porosity  as  about  0.001 . 
The  water  re-entry  pressure  is  given  by  eq  5d 

F^„  =  0.1180/rf^.  (A2) 

For  assumed  values  of  k.  and  b,  eq  Al  yields  the  vein  width  d^.  Eq  A2  then  provides  the 
re-entry  pressure 

We  solve  eq  1 1,  Al  and  A2  simultaneously  by  trial  and  error,  inputting  values  for  k.  and  b 
until  the  observed  time  is  matched  for  a  height  of  z  =  0.043.  This  simplified  scenario  should 
provide  order-of-magnitude  estimates  of  k. 


15 


To  Contents 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302. 
Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid 
0MB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YY) 

September  2001 


2.  REPORT  TYPE 

Technical  Report 


3.  DATES  COVERED  (From  -  To) 


4.  TITLE  AND  SUBTITLE 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


Capillary  Tension  in  Rotting  Ice  Layers 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 


Rachel  E.  Jordan  and  Jeffrey  A.  Stark 


5d.  PROJECT  NUMBER 

61102/AT24 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 

AT24-SP-001 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Engineer  Research  and  Development  Center 
Cold  Regions  Research  and  Engineering  Laboratory 
72  Lyme  Road 
Hanover,  NH  03755-1290 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


ERDC/CRREL  TR-01-13 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Department  of  Defense 
Washington,  D.C.  20314 


10.  SPONSOR  /  MONITOR’S  ACRONYM(S) 


11.  SPONSOR  /  MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

Available  fromNTIS,  Springfield,  Virginia  22161. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

This  study  examines  the  moisture-tension  constituent  relationship  of  ice  near  0°C  by  measuring  its  moisture  characteristic  curve.  In  this 
novel  approach,  we  held  ice  samples  in  a  Tempe  cell  submerged  in  a  cold  bath  and  removed  water  from  the  ice  veins  through  incremental 
increases  in  air  pressure.  Observed  air  entry  and  water  re-entry  pressures  generally  support  an  idealized  flow  path  model  of  alternating  vein 
segments  and  nodes.  Re-entry  of  water  into  the  air- filled  ice  samples  occurred  slowly  over  a  period  of  weeks,  suggesting  low  permeabilities 
on  the  order  of  1  x  10“l^  m^.  We  conclude  that  ice  desaturates  and  rewets  in  a  manner  similar  to  soils,  with  entry  pressures  on  the  order  of 
silts. 
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